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1. INTRODUCTION 



This report documents the results of a six-month study on the pointing 
control of a microwave antenna of the Satellite Power System (SPS). The study, 
performed under NASA contract NAS8-33604, may be considered as the continuation 
of an early effort on the accuracy analysis of pointing control system of 
SPS, which was conducted, under NASA contract NAS8-33065. Primary goals of 
the study encompass three areas, namely, the SPS antenna pointing error sensing 
method, a rigid body pointing control design study, and approaches for modeling 
the flexible body characteristics of the solar collector. 

Accuracy requirements for the antenna pointing control consist of a 
mechanical pointing control accuracy of three (3) arc-minutes and an electronic 
phased array pointing accuracy of three (3) arc-seconds. Results of the 
study, based on the factors considered in current analysis, show that the 
three arc-minute overall pointing control accuracy can be achieved in practice. 

Section 2 of this report deals with the antenna pointing error sensing 
method and its performance. Among the topics discussed are the pointing error 
sensing method, received signals, subarray signal processing, array signal 
processing, effect of noise, and effect of vibrations.' 

Section 3 presents the details of a mechanical pointing control design 
for the antenna. Included in the section are the baseline satellite con- 
figuration and a set of baseline parameters, an estimation of motor rating 
and motor parameters, slip-ring friction characteristics, the effect of 
centrifugal force-induced torque on antenna due to gimbal rotation, the 
effect of variation of antenna moment of inertia due to its changing attitude, 
control system modeling, control system design, system simulation, and the 
effect of noise in sensed pointing error signal. Also included in the section 


is the effect of collector vibration on the accuracy of pointing control. 
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Shortly after the start of the study, It was found that the required 
effort for an analytical modeling of the flexible structure of the solar 
collector was Immense, way beyond the planned effort. As a result, v.'ork 
in this area was redirected into two directions. The first direction is 
to explore the approaches for obtaining an analytic flexible model, and 
the second direction is to estimate the effect of the collector's flexible 
body vibration on the antenna pointing control accuracy. Section 3 is 
devoted to the discussion of the approaches for the flexible collector 
modeling. 

A list of references and two appendices are included for the con- 


venience of readers. 




2. ERROR SENSING METHOD FOR WITENNA POINTING 

In the final report of an earlier effort [1] a method was proposed 
for developing angle error signals to be used by angle servomechanisms 
to point the antenna toward the ground station. In this effort, the 
proposed method has been shown to be feasible. In addition, a normalization 
technique, similar to that used in angle tracking radar, has been added 
that produces error signals that are independent of received signal magnitude. 
Finally, sane results related to affects of vibration and noise have been 
derived. 

2.1 Sensing Method 

The basic concept is illustrated in Figure 1. Two intersecting linear 
arrays formed from subarrays of the main array are used to develop error 
signals. The array normal to axis 1 is used to develop an error signal 
proportional to the error that the antenna pointing direction makes relative 
to the desired direction (toward the earth site) in a plane for which axis 

1 is normal (error about axis 1). The other array generates a signal 

proportional to the (space) orthogonal angle in the plane for which axis 

2 is normal. The two error signals are used by angle servomechanisms to 

turn the antenna about axes 1 and 2 to cause the errors to approach zero. 

Each array is presumed to have 2N + 1 * M subarrays with the conmon 

subarray at the antenna's center serving as the "reference" for the two 
arrays. In general, the two arrays do not have to have equal numbers of 
subarrays; however, since accuracy requirements are likely to be the same 
in both coordinates, there seems to be no reason to not make the two 
identical, as is assumed in following discussions. 
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Subarrays applicable 
to coordinate 1 



Figure 1. Subarrays applicable to ginbal servo error generation. 


2.2 Received Signals 
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We define antenna-centered coordinates as shown In Figure 2 where 
X and y are axes In the plane of the sntenna's face. Axis z Is the 
antenna's true pointing d"!rect1on. A wave arriving from the earth site 
makes an angle 6 from the : (boresight axis) direction; It has an 
"orientation" defined by an angle $ in the xy plane. Using these definitions 
we find the signal received by a typical subarray located along the x 
direction. Only the x axis array Is analyzed since the other behaves in 
a similar manner. 

The error angles sensed by the system in Figure 2 are 6^ (In the xz 
plane) and (In the yz plane). The various angles are related by 


tan(6j^) ■ tan(e) cos($) 

(1) 

tan(e ) ■ tan(e) s1n(o) 

(2) 

from geometry, 

A typical subarray pattern in the SPS system will have a 
illumination function, have a rectangular shape (side-lengths 
a one-way voltage pattern, denoted G(e.^) by 

uniform 
X and Y) and 

G(e,(^) ■ Sa [“• sin(e)cos(<>)3SaC^ sin(e)sin(^)] 

(3) 

where 

s ,(0 . 

(4) 

and X is wavelenght. If we now recognize that 'ei<<l for any reasonable 

SPS system, then |6 |<<1 and |<<1, so (1) through (3) give 

X y 

^ e cos(o) 

(5) 

e > e sin c 

V ^ 

(6) 
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Figure 2. Antenna-centered coordinates. 
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6(ei») % S«[(’tX/^)e cos(t)] SaC(^Y/x)s sln(f)] 

• SaC^rXe /x]Sa[itYe /a] - G(9 .s J (7) 

A y A / 

Now considtr the array of subarrays along the x axis as shown in 
Figure 3. Subarray locations along x are nX, -N<n<N, and they have dis- 
placement in the z direction of z^ and tilt angles (tilt in the other 
coordinate Is ''n voltage produced by the nth sub- 

array by a wave having error angles e and 6 , then 

where is the voltage produced when the subarray points toward the 
pilot source without tilt error and Is the phase of the arriving 

pilot signal at the reference (n«0) subarray. Since l6j^l<<l. 

2TrnXe 2ftz 

2.3 Subarray Signal Processing 

We next show that desirable error sensing occurs if the signal pro- 
cessing of Figure 4(a) is available at each subarray in the line array. 
Straightforward analysis gives 


e_^*^I A. 2rrnXe 2*(z -z-) 

«n 

2^nX6 2^(z -Zfx) 

j;.. 


( 10 ) 


^intcriOfJ 

/ 

i 


/ 



X • n 


Figure 3. Subarray locations. 
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The signals and are assuned conveyed to a central processor for 

the array, as shown in Figure 4(b). 

2.4 Array Signal Processing - No Errors 

Suppose no tilt or displacement errors are present. Then 

9 „=0, and 2„*0 for all n. Now (10) becomes 
yn n 


e A 

" > ■ (coA[n(¥ 

( ) 4 X y A. X q 

®n 

+cos[n(~ 9x"®q^^^- 

The signals e^'*’^ and e^'^ in Figure 4(b) become 



Aq sin[^(^ 


sinC^-C 




where 




} 


M*2N+1 


and we have used the known series 


N 


. p2N+1 

sin[-^ 


[ cos(n6) ^0 

n="N sin( Z^) 



( 11 ) 


( 12 ) 


(13) 


( 14 ) 


Next, define a subarray beamwidth (in rad) between -392dB points 


n 


as Sjytj* Siven by 


and 




®sub 


9q/27T. 


Then (12 becomes 


(15) 

(16) 


e<*>, A„ si„[Sl_(e +e )] 




sin[^ (9 +e„)] 


•0 


sub 


X q 


s1n[^(9 -? )] 
^ ^sub ^ ^ 

' sin[| — (e -9„)] 


'sub 


X q' 


(17) 


Now if we define an "array pattern" G (e ) by 

o X 


sin[ 


Mtt 




M fr 9 

ej sin(^ ^) 


sub " °sub 

7T e ^ TT 0 

/ x> 


for !e,l«6s„b/r 


M sin(^) 
®sub 


0 


sub 


(18) 


then 


e^^^ AM 

(_)^*“r ‘^(®x* ®y^ ^'^a^V®q^ - '^a^V®q^^- 


( 19 ) 


A sketch of G (e^) is shown in Figure 5 for M=3,5, and 9. The function 
is symmetric about 2=0 and + ft and about z-+^l2. Thus, there is another 
peak in the function at 2 =+ft. For M large G.(e ) approaches a sin(x)/x 

* o X 
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function for small magnitudes of x. The responses (peaks) at never 
show up in the response (17), however, because G(e , 6 )«0 for z»’tX 6 /x • 

A A 

ire /6 . ■ 'T from (7). 

X sub ' ' 

Next, define array sun and difference patterns by 


V«x> - 

respectively. Equations of (19) become 


I e^^^ « Aq M 6(e^, 6y) G,(e^) 

I e^'^ » A- M G(e , 0 ) GJe ). I 

I 0 ' X y' A x i 

' I 

The response e^ in Figure 4(b) becomes 

! 


( 21 ) 


( 22 ) 

(23) 


(24) 


The behaviors of the sum, difference, and "normalized difference 
pattern" G. (\)/Gr(e„) are shown in Figures 6-9. In these figures we 
define an "array basic beamwidth" by 

®B " Sub^^ ‘ 

Figure 6 shows that the sum pattern can be of the usual single-peak form 
or double-peaked. The transition takes place when s^/Sg ■ 0.6625. Figure 
7 shows typical difference channel responses with Figure 8 showing that the 
slope of the response at the origin is maximum v/hen ?q/6g is about 0.66. 
Figure 8 also shows the loss in gain of the sum channel as a function of 
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Sq/9g. Finally, Figure 9 shows the normalized difference pattern. From 

these data, which apply for M larger than about 9, It Is clear that ?„/3o 

q ts 

near 0.5 should be a good system choice based on small sum channel gain 
loss and beamwidth broadening, large error slope, and linearity of the 
normalized difference pattern. 

For 9q/9g *0.5 the slope of the difference pattern on boresight Is 


slope 


“3“ e. 


ex-0 


»4 , -4 M 

" ®sub 


(26) 


The slope of the normalized difference pattern on boresight Is 


slope 



r^i 


e. -0 



-2 M 
®sub ' 


(27) 


2.5 Effects of Molse 

Assume each subarray output contains a bandpass gausslan noise N^(t). 
Such noises can be described by [2] 

• N,„(t) cos[«p,,^t t] - N„„(t) t] (2S) 

where fljp(t) and N^^(t) are Independent, gausslan random variables for any 

t and a given value of n. They are also Independent between subarrays 

(different values of n). Furthermore E[N^^(t)] ■ E[N^^^(t)] » E[N^^^(t)] * 

Is the noise power (assumed the same at each subarray). 

Straightforward analysis of Figure 4(a), assuming the LPF's pass 

the full noise bandwidth, gives the noises present in and e^”\ 

n n 

denoted and respecti vely. 



[N^^(t) COS(9j) + N 


on 


(t) 


sin(e-.)] cos'ne_) 

' 'A 


( 29 ) 
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• [Nj„(t) Sin(e^) - cos(e^)] sln(ne ) 


(30) 


Here ^pi‘io^(t) ■ t ■*■ e-f. where 8j is an arbitrary pilot source 

phase angle. 

After sirmlng the various noises as Implied In Figure 4(b) the 
noises on the signals e^"^ and denoted N^“^(t) and N^”*’^(t), respectively, 
are found. The powers In these noises follow taking the mean-squared 
values: 

4 “ - 0^ c| + 7sin(ej'^ 3 (31) 


E[(N(-)(t)}2] . 


(32) 


Figure 10 plots (31) and (32) versus ^g/6g ■ M 8^/2^. Clearly, for squint 

— 2 2 

angles of Interest (e^ ■ ®Z " output noise powers 

are equal ; 


2.6 Effect of Vibration (Displacement) Errors 

We assume that there Is symmetry in the vibration of the overall 
array so that the geometry of Figure 11 applies. Since subarray positions 
are 

■ n X, -N < n < N, (34) 

then 


2n ■ z(Xj^) ■ A2 cos[(2K-l) IT X n/D]. 


(35) 


For small displacements lz « D, as Is the case in practice, it is readily 
shown that the vibration induced tilt is 

5 . 


sin[(2K-l) tt X n/D]. 


(36) 



0 


X 


X 


^(t) m 6^ l(2K^t)Tr)f/V] 

D*Aa/t£a/^a V/Aj^ierM/t 
K * VtBkAnoN tAooe, • 


J»* * ^5 c^.f^air-oi^jr h/p 3 

5*^ l(iK-0>X *»/^J 

X„ i V X 


Figure 'i. Displacement on the antenna over a line 
through its center (origin). 
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NumerlcaViy, If • 1 m, 0 • 1010 m, X ■ 11. o4 m, and .x • 0.1224 m, 
then ^ 0.89 so tilt errors may be nearly as large as the 
subarray's beamwidth. 

We evaluate vibration In the x direction when • 0, ■ 0, and 

when 

iiz ■ dZ cos(wyt) (37) 

where -iy/2*r Is the frequency of the vibrations. From (10), using (7), 
the signals e^"*"^ and e^‘^ In Figure 4(b) are 


e 

e 


(*) 

(-) 


A. N 

-w- I Sa[TrX(e +6 p)/>,] 
^ n>-N 


2irXe 2ir(i -Zq) 

{cos[n(— 2-^ ej * -'^-^3 


2ttX6jj 

+ cos[n(— 6^) * 


A 


) 


(38) 


For relatively small values of the system's response will be 

However, the true response without vibration will very 
nearly equal the s1;ioe when ■ 0 times *9^^, so that the "error" in radians 
produced by vibration, denoted 9y, will be approximately 



after (27) Is used. If normalized to the jrray basic beamwidth 6g we have 


9 

e 


V 

B 


a 



( 40 ) 


while (38) becomes (for « 0.5 Sg) 
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•{cos[n( 
+ cos[n( 


Zv 9 

wt; 

Zir 6 

WT^ 


X 


X 




X ■* 
X 


(41) 


An examination of (41) shows that * 0 if * 0. Thus, for 
no antenna pointing error (e^ = 0) the presence of vibration produces no 
error , reqardles : £* the vibration mode , so long as it i^ symmetri c . 

To study the error 9y/6g when ^ 0, (40) was computed on a digital 
computer for various, values of M, K, and when peak vibration dis- 
placement in each case is 1.0 m. Figure 12 illustrates curves of 6y. as 
a percentage of array basic beamwidth 9g versus the number M of subarrays 
in the array. Curves are shown for pointing error magnitudes |e /ed of 
0.05 and 0.1 for two vibration modes (K=1,2). The curves assume X=11.6A m, 
0=1012.63 m, and x=0.1224 m. It is clear that sane combi nationsof vibration 
mode and pointing error give rise to large errors due to vibration. 

For example, with M=7 subarrays in the linear array and mode 2 vib- 
ration (K=2), when pointing error 9^ is 5 % of the array's basic beamwidth 
9g, the maximum error due to vibration of peak amplitude 1.0 m is 7.32% of 
9g. Figure 13 plots 6g versus M and gives 9g = 5.16min for M=7. Thus, 
maximum error due to vibration is 0.0732 (5.16) = 0.38min. Since the 


pointing error is 0.05 (5.16) = 0.268min, we find the vibration error even 
exceeds the pointing error to be corrected by the servo system. The vib- 
ration error cycles with time, and, for the example taken here, is illustrated 
in Figure 14. 





i 




m 


n 
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Finally, we note that vibration will have an affect on signal -to-noise 
ratio of the composite signals and e^”^ The affect on e^'*’^ is 
of great concern since signal -to-noise ratio for this signal should be 
maintained above 10 to 100 for good performance of practical normalizing 
circuits. Figure 15 shows the maximum loss in this signal -to-noise ratio 
for various M, K, and 9j(/e0 when az » 1.0 m. 




3. RIGID BODY ANTENNA POINTING CONTROL 


In this section, an antenna pointing control system is modeled and 
then designed assuming that both the solar collector and antenna are rigid 
bodies. The purpose of this design is two-fold. First, the design will 
generate technical data such as motor requirements, system response time, 
effect of noise and system parameter changes, and the stiffness of the 
control loop. Second, by examining time constant of the control system 
and the mode frequencies of the flexible collector, the effect of the 
collector vibration mode on antenna pointing accuracy can be estimated. 

The goal of the mechanical pointing control is to maintain a mechanical 
pointing accuracy of three arc-minutes. 

3.1 The Baseline Satellite Configuration 

Figure 16 depicts the baseline satellite configuration. The dimension 
of the solar energy collector and the position of the microwave antenna 
are shown in Figure 16(a). The size and the general shape of the yoke- 
antenna combination is shown in Figure 16(b). Figure 16(c) is included 
to show the size of the rotary joint connecting the yoke to the solar 
collector. 

The mechanical control of the antenna pointing is done by rotational 
motions about two axes. One rotational motion is about the Zy axis which 
is along the longitudinal center line of the collector as shown in Figure 
16(a). A second rotational motion is about the X„ axis which is along 
the direction of hinges joining the antenna and the yoke as shown in 
Figure 16(b). 

Two coordinate frames are used in the sequel, one is the antenna 
frame (X^,Y^,Z^) and the other is the yoke frame (Xy,Yv,2y). 

29 



16. Baseline configuration for SPS satellite. 
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The physical parameters of the satellite are collected and listed 
here. For the collector: 

Mass: ■ 16x10^ kg 

Dimension ■ i ■ 17,000 m 
a ■ 4,200 m 
b ■ 600 m 

moments of inertia ■ 

(Ic)z^ - ^ (a^V) « 19.5x10^2 kg-m^ 

dc)j(^ * j - 1.54x10^^ kg-m^ 

For the combination of antenna and yoke: 

Mass * ■ 13x10® kg 

Diameter: D * 1 ,100 m 

Moments of inertia: 

(Ia)z = » ^MD^ » .983x10^^ kg-m^ 

“ •g- MR^ « .983x10^^ kgV 


Other physical data include: 


slip-ring diameter* 350 m 
motor shaft diameter » .1 m 



The gear-train is used to connect the motor shaft to the rotary joint, 
the rotary joint has a diameter of 350 m and the motor diameter is .1 
a gear- train ratio of 10"^ is reasonable. 


Since 
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3.2 Slip-ring Friction Characteristics 

Conduction of electric power from the solar collector to the antenna 
yoke and from the yoke to the antenna Is done by using slip-rings. The 
friction torque at each rotary Joint Is mostly due to the slip-ring friction 
torque which is a nonlinear function of speed of rotation. The nonlinear 
slip-ring friction characteristics, as provided by NASA/MSFC, are shown 
In Figure 17. In the figure e represents the relative angular position 
between two sides of the rotary joint. The arrowheads Indicate directions 
of motion. Notice that the friction torque has the saturation values of 
+ T_.„ » + 10® N-M. In the region between the saturation values, the 
friction troque depends on e in a way resembling the compliance torque 
The compliance is the slope Kp which Is -10® N-M per degree. The nonlinear 
friction torque exhibits the nature of shifting hysteresis. 

Figure 18 gives two analog electronic circuits which r ,n be used to 
simulate the above friction characteristic' on an analog computer. For 
system simulation using a digital computer, however, neither one of these 
circuits are convenient to use. To provide a convenient way for later 
digital simulation, a digital simulation flow-diagram is developed as 
shown in Figure 19. 

3.3 Estimation of Motor Parameters 

It is assumed that d-c motors are used as torques for antenna pointing 
control. The rating fcr the d-c motors are estimated here. Assume a maximum 
perturbation angle for the antenna to be [3] 

l^el * 1 degree * .0175 radian 

nifl X 

The moment of inertia of the yoke-antenna load is 
I^ » .983x10^^ kg-m^ 



friction characteristics exhibits shifting hysteresis 




n 









t 

OUTPUT Tic 

Tk-i 


Figure 19. Digital simulation for slip-ring friction 
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Giver a aear-traln ratio of the reflected maximum angle perturbation 

and the moment of Inertia of the load at the motor shaft Is 

• 9830 kg-m^ 


1^6*1 ■ — lAfli ■ 17S radians 

' 'max a ' ‘max ' ra«ian> 


(see discussion In Scctlm 3.6). 

The maximum slip-ring friction torpue 1$ 


Itplmax • N-M at :he ring 


■ 0x10® • 100 N-M at the motor shaft 


The Initial motor torque for a maximum step perturbation Is 


I'ol • ‘a i« 

max 


U l^®'Lax _ 1.72x10® .. 

n ■ ■ ■ ' - m 


where T Is the system time constant. Therefore the total maximum torque is 


max * ‘"'Fornax * '’0‘max 


100 + 


1.72x10' 


An estimate of motor power P Is now given 


P • Average of [A6{t) T(t)] 


< U ' 1 It! 

- ' ‘max ‘‘max 


jaS' 

' ‘max ' ‘max 




1.75x10^ T^+3x10® 


watts 


r. 


Ul 
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23.46 T^+4.02^x^0^ 


Under the nominal condition the steady-state motor power is to overcome 
the slip-ring friction torque. This power is 

A list of P a:n ve 'sus T is shown in Table 1. A baseline rating is 
chosen to be: 

P « 1681 watts = 2.253 HP 
T = 578 N-M 
T * 60 seconds 

Table 1. Motor Ratings 


T (sec) 

T (N-M) 

P (Watts) 

10 

17293 

301750 

30 

2010 

11694 

60 

578 

1681 

100 

272 

475 

300 

119 

69 

600 

105 

31 

1200 

101 

15 

1800 

100 

10 


Adopting the baseline motor rating, motor parameters are now estimated. 
Take the effective d-c voltage to be 

Vdc = ^20 V. 


The average motor armature current is then 




r 


L. 


V 


dc 


1681 

120 


14 A 


The motor torque constant is 

578 


't I. 14 A 


Using the MKS system of units, the motor emf constant equals the motor 
torque constant in value, giving 


Kg * 41.3 V-sec 


R » ■ .43!! 

V 


The result is summarized in Table 2 for the ease of future reference. 


Table 2. Motor Parameters 

Rated Power 

1681 W 

Rated Voltage 

120 V 

Torque constant, ICj. 

41 3 Ml 

Emf constant. Kg 

41.3 V-sec 

Armature resistance, R, 

a 

.43 n 


3.4 Effect of Centrifugal Force Induced Torque 

Under the nominal condition the yoke-antenna assembly is rotating 
about the Zy axis of the SPS satellite at nearly constant speed. The 
rotation generates on each mass element of the antenna a centrifugal 
force which is directed perpendicularly to and away from the Zy axis. 
This force produces a torque about the axis of the antenna, which may 
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Assuming a 5« motor power loss, the motor armature resistance is computed 
by 


39 

be viewed as a load disturbance torque for the pointing control system. 

This torque is analyzed with the help of Figure 20. 

Lat aM be the increnental strip of antenna mass as shown in Figure 
20. Let be the nominal angular velocity of the yoke-antenna assembly 
about the Zy axis, which is very close to the Earth's rate of rotation. 

The incremental torque about the axis due to the centrifugal force on 
Af-1 is 

= {AM)(acceleration)(torque arm) 

= (aX^ Mp)[2l sins)(2^ cose) (46) 

where o is the area mass density of the antenna, which is assumed constant. 
Integrating over the entire antenna plane gives centirfugal torque as 

f ^ 2 

Tc = j ^ J ^ \ 

area 0 

2 ^ ? 2 

= 2 o sin2e | /R^-2^ dZ^ 

'o 

where R is the radius of the antenna. After the integration the total 
centrifugal torque is obtained as 

Tc = ^ 0 sin2e (47) 

Maximum torque occurs when e*45°. The Earth's rotation rate is 

-5 

= 13 deg/hour = 7,27x10 rad/sec 

Using the baseline parameter values, R = 550 m and 

3 = = 13.68 kg/m^ 

^ 77x550"^ 



* 


Centrifugal force Induced torque. 




Therefore the maximum torque is 


‘c ‘max 


= 2598.2 N-m 


(48) 


Comparing this torque to the maximum slip-ring friction torque, 
which is loS-m, the former is less than .3« of the latter. The effect 
of the centrifugal force induced torque therefore needs no special 
attention, but may be lumped into the slip-ring friction torque as a 
perturbation. It is expected that such a small perturbation will not 
cause problems for the antenna pointing control. Computer simulation 
will be used to verify tnis fact. 


3.5 Effect of Variation in Antenna's Moment of Inerita 

Since the antenna disc may rotate about the axis as indicated in 
Figure 21, its angular position e about this axis may change. As a result, 
the antenna’s moment of inertia about the Zy axis may change. It is 
desired that such change will not have adverse effects on the accuracy 
of the pointing control. 

Let the moment of inertia matrix of the antenna expressed along the 
antenna frame (X^,Y^,Z^) be 


I 


X 


0 


0 


I 


0 


0 


I 


y 


0 


0 


I 


2 


Referring to Figure 21, the antenna frame and the yoke frame are related 
by an angle e. The antenna's moment of inertia expressed along the yoke 
frame (Xy,Yy,ZY) is 
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“T 

“ 


— 

1 

0 

0 


I 

X 

0 

0 
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0 

0 

0 

ce 

se 


0 

ly 

0 


0 

ce 

-se 

0 

-se 

ce 


0 

0 



0 

se 

ce 

_ 












0 ce^+I^ se^ 

0 -ly cese+I^ces8 


0 

-I^cose+I^sece 


(49) 


where c0 = cose and se = sine. In particular, the moment of inertia about 
the Zy axis is 

i; = ly se^I^ ce^ (50) 


Note 



1 2 

M^R = max, ant. mom. of inertia 
1 2 

^ M^R = min. ant. mom. of inertia 


Taking the derivative of (50) with respect to e 

^ = 2 se c9 - 2 se ce = (ly-I^sin 2e 

Therefore the variation of the moment of inertia normalized with respect 
to 1 2 is 

-j — = (j - 1) sin 2e A9 = sin 2e ie (51) 

z z 

For e=45®, (51) is at its maximum. For a 9=1 degree = ,0175 radian, the 
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specified attitude control accuracy for the SPS satellite, 


r 

I 


.0175 « 1.75% 


(52) 


max 


A plant parameter variation shown in (52) is too small to cause adverse 
effects in pointing control, and can be ignored in the design of pointing 
control. However, computer simulation will be performed to confirm this 
observation. 


3.6 Pointing Control System Modeling 

An all -attitude pointing control capability of the SPS antenna is made 
possible by allowing two degrees of rotational freedom for the antenna disc 
as shown in Figure 22. The first rotary joint allows the rotation of the 
antenna-yoke assembly with respect to the collector about the Zy axis, while 
the second rotary joint gives rotational freedom for the antenna with 
respect to the yoke about the axis. It is assumed that the coupling of 
the antenna's motion about these two axis is negligible. This is due to 
the fact that the motion of the antenna is very slow and that the angular 
p . "turbations about the two axes are very small. Therefore the dynamics of 
the pointing control system for each rotary joint may be separately modeled. 

The First Rotary Joint Let and be the absolute angular positions 
of the antenna-yoke assembly and the collector, respectively, about the 2y 
axis. Their dynamics are given by 


‘a ' - ’F 




(53) 


( 54 ) 
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where and are moments of inertia 
the collector, respectively; t' is the 
motoi"; and is the frictional torque 
compliance and other kinds of friction 


of the antenna-yoke assembly and 
control torque originated from the 
of the slip-ring. Rotational 
are assumed negligible. Adding 


(53) and (54) gives 


A A c c 


(55) 


Define the ratio of moments of inertia 


6 


*A 


(56) 


Then (55) gives 


ee, 


(57) 


The relationship among the generated torque t, applied voltage V, and 
the shaft speed se for a separately excited d-c motor is 

V-Kpse 

T * IC^ — (58) 

where K_, K- and R are defined in Section 3.3. Noting the presence of 
I b d 

a gear-train with a ratio a, the following relationships exist. 




(59) 


T 


a T 


(60) 


Ignore the nonlinear friction torque Tp in (53) for the time being. Then 

I, (61) 

Eliminating 9 ,t and t' from (53), (59), (60) and (61), yield 



[(»^I,)S^ 






a 


(62) 
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Let 



(63) 


which represents the reflected moment of Inertia of the antenna-yoke 
assembly seen at the motor shaft. The transfer function of the first 
rotary system Including the motor Is therefore 


T 


m 


a ICp 


S[S+T^- 


-( 1 + 6 )] 


(64) 


For later convenience, define 



(65) 


which Is the reflected absolute attitute of the antenna-yoke as seen by 
the motor. Using (64) In (65) gives 


*P ‘ T “ 


Me 

s[s+yp[~ (1+6)] 


( 66 ) 


Eq. (66) Is the transfer function of the plant at the first rotary joint. 
Notice that the effect of the reaction motion of the finite mass collector 
Is represented by 6 In the equation. For an Infinite mass collector, 6*0. 
A block diagram for this plant is given in Figure 23(a). 

In sections 3.1 and 3.3, values of various parameters have been 
obtained as: 


* 9830 kg-m^ 
K.J. « 41.3 

* 41.3 V-s 


.43 n 



a 
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Quantity s as defined by (55) has a value of 

e . . .06 

19.5x10'^ 

M1th numerical values Inserted, (66) becune' 

" T ’ s”ts+'.4z4) 

In this model the nonlinear friction was not Included. The effect of 

the nonlinear friction will be treated as a torque disturbance In the control 

systen design. 

The Second Rotary Joint For the control about the axis at the 
second rotary joint, the collector's moment of Inertia may be assumed 
infinite. This Is demonstrated as follows. Referring to Figure 24 and 
using the baseline data given in Section 3.1, moments of Inertia about 
axis are: 

For the antenna: 

”a>x. ■ I 

For the collector 

• 394.3 

Thus, a rotation of antenna about by 1® causes a rotation of collector 
about X^ less than 3 arc-seconds. For a 3 arc-m<nute mechanical pointing 
accuracy, the assumption of infinite moment of inertia contributes only 
.11 arc-second error. 


4 # 



Moment of inertia 
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The following are further assumed; (1) the same type of d-c motor 
will be used at the second rotary joint as the type used at the first 
rotary joint. (2) Same gear-ratio. (3) Mass of yoke is small as compared 
to that of antenna. Under this condition the transfer function for pointing 
control at the second rotary joint can be obtained from (66) by setting 6*0, 
yei Iding 




"A 

T 


TT-T T ' 


:[s+y 




. R- ^ 
A a 


( 68 ) 


where is the reflected antenna angle about the axis as seen by 
the motor. The actual antenna angle is 

?A * 

With numerical values inserted, (68) gives 

r _^A . .01 

S ■ T ■ sTsOO^ 

A block diagram of the plant representing the second rotary axis is given 
in Figure 23(b). 

3.7 Pointing Control System Design 
Design Background 

Under the assumption of negligible dynamical coupling between 
the two axes of rotation, the two associated control channels can be 
separately designed. The need of pointing control is to counter-act the 
effect of external disturbances and changes in plant parameters. External 
disturbances include gravity gradient torque, solar pressure torque, and 
collector notion coupled through the rotary joint. Change in plant para- 
meters include the change in the antenna's moment of inertia due to change 


(69) 


(70) 
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in antenna attitude, change in motor armature ressistance caused by thermal 
conditions, aging effects, change of structural shapes and dimensions 
due to temperature changes, and other-'-.. Most of the disturbances are con- 
stant or very low frequency cyclic. Parameter changes are also very low 
frequency in nature. It is assumed that the frequency range for distur- 
bances and parameter changes is below .001 hertz (15 minutes per cycle). 
Effects of instrument noise will be ignored because they will be well 
attenuated by the low-pass characteristics of the system. 

Under the nominal condition, the antenna-yoke is rotating at the 
first rotary joint with a constant rate of one revolution per 24 hours. 
Magnitudes of disturbances at this joint are assumed small so the direction 
of rotation is not reversed. The second rotary joint is stationary under 
the nominal condition. With disturbances added, it may rotate in either 
direction. 

Comparing (67) and (70) shows that they differ only in time constant 
of an amount of less than 5 %. Many well designed control loop for one plant 
should therefore perform well for the other. In view of this, a control 
system design will be made for (67), the plant for the first rotary joint. 

Figure 25 depicts an antenna pointing control loop for the first rotary 
joint. An ideal gear-train is assumed. The pointing error is sensed 
by the microwave error sensing method presented in Section 2.1. An 
algorithm is used to generate, from the error signal, a command signal 
for pointing control. Three ways of generating command signal have been 
developed previously. [1] For the convenience of design the transfer 
function between error e and command signal may be taken to be unity. 

The present control system design is to determine the block labeled "network 
and amplifier gain" . 






Antenna pointing control 







54 


The Design Except for the 3 arc-minute mechanical pointing accuracy 
requirement, no other system performance specifications were given. 

Linear quadratic optimal control design will be adopted. Since the plant, 
represented by (67), has single-input and single outout, frequency domain 
optimization by way of the root-square-locus technique will be used. [4] 
Because of the constant rate rotation under the nominal conditions, the 
control system for the first rotary joint should be a follow-up system. 
Follow-up design will also be used for the second rotary joint to provide 
flexibility for the possible need of follow-up operations. 

Begin with (67), the plant transfer function* 

r - -01 

'’p ■ s(s+.424) 

which has poles at p.j=0, P 2 * 0 , and p^=-.424. Open-loop poles for the 
root-square locus plot are 



There are no open-loop finite zeros. The root-square loci are shown in 
Figure 26. Choose the closed-loop root-square poles at 


Then 


Q., S Qj * .8 + j 1.48 
Q 3 « -.18-0^-02 = - 1.78 

-q.| & «= --S + j 1.48 

.Q 3 = 1.78 


ir 

An integrator 


is added to fiiake the 


system type-2. 
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The closed loop poles for the control system are given by 


« - /-Q^ » 1.297 1 120.8 


Compute 


. . /-Q^ - 1.297 1 -120.8 


^3 ' - >^-03 * -1*33 


q.jq 2 “ 1.682 

q.|+q 2 » -1.328 

-L/ - - .79 
qi92 


The reciprocal of the closed loop zero is 

<’1 ^2 '’a 


T= . (J-.-L.JL) = . (!l!2i.J-) 


» .79 + .75 * 1.54 


The closed-loop transfer function is 


(1+Ts) 


fUl- f) 


‘'2 


1+Ts 


(T . J_)(i . iL-li s + — ^ s^) 

^1 ^ 1*^2 


1 + 1.54 s 


(U.75s)(H.79s + .595s^) 


1 + 1.54s 


T 


1 + 1.54s + 1 . 18s“ + . 45 s 


(71) 


(72) 
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Step response simulation was performed to test (72). The overshoot 
was found to be more than 40^: not a good response. Repeated trial - 
and-error adjustment led to the following closed-loop transfer function. 


H » X ' 2 3 

1 + 3.164s + 3.082s'^ + .595s'’ 

This system had a step responce overshoot of 18.7* and a time-to-peak 
of 3.4 seconds. The corresponding open-loop transfer function Gq was 
obtained from 


giving 



.3245(1 3.164s) 

s^(l + .193s) 


( 73 ) 


(74) 


(75) 


Two compensation configurations were used as shown in Figure 27; 
the first is a series compensation and the second if a series compensation 
with rate feedback. The second configuration offers better counter-action 
for load disturbance and plant parameter variation. 

For series compensation: 

. _ ,, n * 3,164s)(l ♦ 2.36s) 

Gc ■ g-- 13.75 J s(l * .193s) ' 

For series compensation with rate feedback: 


choose K.^ » 475.76 * 476 
then G, . 


(77) 

(78) 


A block diagram for the simulation of the designed closed-loop 
system is given in Figure 28. The dashed line part is used only for 
series compensation with rate feedback. Figure 29 shows the step response 
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of the system under the nominal conditions. The reference Input Is 175 
radians for motor, corresponding to 1* for the antenha. It should be 
mentioned that the system Is not expected to be subjected to a step 
function type of cotmiand. Step response simulation Is used merely to 
reveal the system's transient characteristics. 

System Sensitivity The response sensitivity of the designed system 
was studied by way of computer simulation. Results are given in Table 
3. As a whole, the response did not change much for parameter change 
as high as 20%. Also, the superiority of the series compensation with 
rate feedback is evident. 

The Constant Drag Torque Since the .first rotary joint does not change 
its direction of rotation the nonlinear slip-ring friction torque appears 
as a constant drag torque with a magnitude of entering the system 
at the point as shown in Figure 30. Results of a step response simulation 
are given in Table 4. One sees that the control is capable of nulling 
the effect of constant drag torque very effectively. A CSMP listing of the 
simulation is attached as Appendix A. 


Table 4. Response of System with Prig Friciton 


Overshoot n 

Time to peak tp 

Error at 4 t: motor angle 

P 

antenna angle 


21 % 

3.3 sec 
.13 rad 
.04 arc -min 


Theoretically, the steady-state pointing error due to a constant 

dr?.g torque should be zero. From Figure 30, the response of the system 

A ■ 

jutput due to a constant disturbance torque “ ■j ’s 


Table. 2. System Sensitivity. 


Series 

compensation 


Series 

compensation 
with rate 
feedback 


System 

condition 


Nominal 


Kj/R 10 S up 
• « 


’’ 10 % down 


" 20 % up 


" 20 % down 


Nominal 


" 10 up 


" 10 % down 


" 20 % up 


Overshoot 

n 

Time 
to peak 

*p 

Settling 

time 

19 % 

3.4 $ 

9.0 s 

16 % 

3.2 s 

9.6 s 

■■ 

3.6 s 

O 

GO 

CO 

15 % 

3.2 s 

9.8 s 

24 % 

3.8 s 

8.6 s 

ma 

3.2 s 

8.8 s 

19 

3.2 s 

8.8 s 

20 % 

3.2 s 

8.8 s 

1C % 

3.2 s 

8.8 s 

20 % 

3.2 s 

8.6 s 













































constant drag torque 
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s;($) 


2 S 

s + [G^(s) + Kg(l + 3)s] 


The steady-state value of ls 


(79) 


6^ (t—) ■ se^(s-o) ■ 0 (80) 

where the property G^(s-»0) • - has been used. 

The Nonlinear Friction Torque The second rotary joint has little 
rotation under the nominal conditions. To counter anamalous effects, 
it may rotate in either direction. Thus the shifting hysteresis type of 
nonlinear friction torque, as portrayed in Figure 17, may come into play. 

This friction torque is a nonlinear function of the speed of the motor, 
and it enters the system in a way as shown in Figure 31. The nonlinear 
function is designateo 6p. can be simulated according to the flow 
diagram of Figure 19. Mote that all load quantities have been reflected 
into quantities as seen by the motor; the value for Kp in Figure 19 should 
be handled similarly. Using the data in -igure 17, 


•V- 


10^ antenna 

Deg 


57.3x10' 


N-M 


10 


rad 


- .573 at motor 


max 


at antenna 


10° 

» TOO N-M 

10^ 


m 


at motor 


Nonlinear friction torque 



ion torque 
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To reveal the effect of nonlinear friction torque, a square-wave may be 
used as reference input for system simulation since rate reversal is 
needed to exhibit the shifting hysteresis characteristic of the nonlinearity. 
In order to drive the nonlinear friction torque to its saturation level 
the amplitude of the square-wave should be greater than 1® in terms of 
the actual antenna angle. An amplitude of 2® is chosen. The corresponding 
amplitude for the motor angle is (2xl0^)/57,3 * 349.04 radians. A simulation 
was performed in the manner described. . Results of simulation are listed in 
Table 5, which shows the adequancy of the control. Appendix B contains 
the CSMP listing of the simulation. 


Table 5. Square-Wave Response of the System having 
Shifting Hysteresis Type of Friction Torque 


Overshoot n 
Time to peak tp 
Error at 4 tp: Motor angle 
Antenna angle 


10 % 

3.19 sec 
.14 rad 
.05 arc-min 


The Centrifugal Force Induced Torque Recall from the discussion 
in Section 3.4 that the maximum torque on antenna caused by the centrifugal 
force is less than .3% of t which is 100 N-M as seen by the motor. Insert 
T = 1 N-M, which is 1% of T . , into the last simulation program and repeat 

w fHoX 

the simulation. The results are nearly the same as those in Table 5, 
showing that the pointing control is capable of nulling the effect of 
satisfactorily. 
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3 -8 Effect of Noise In Sensed Pointing Error Signal 

It has been shown that, based on the set of given assumptions, the 
designed mechanical pointing control system more than adequately meets 
the 3-arc -minute requirement for pointing accuracy. If the microwave 
sensed pointing error is noise contaminated, however, the situation 
could be quite different. As discussed in Section 2, the noise in sensed 
error is low frequency in nature, which can be as low as d-c. Under tnis 
condition the noise directly dictates the accuracy of the overall pointing 
control system. At this moment, it is not practical to arrive at a 
numerical value for the noise level. The remark here serves to emphasize 
the possible adverse effect of the noise. 

3.9 Effect of Collector Vibrations 

It has been discussed in Section 3.7 that the frequencies of vibration 
for the collector is less than .001 hertz which is considerably lower than 
.2 hertz, the approximate bandwidth of the pointing control system. Therefore 
one expects errors caused by the collector vibration be effectively 
diminished by the pointing control. To support this claim, computer sim- 
ulations were performed as follows. 

Pointing error caused by collector vibration is equivalent to the 
error of the closed-loop control system excited by a sinusoidal reference 
input. Reference inputs having frequencies of .001 hertz (about 16.6 
minute/cycle) and .01 hertz (about 1.66 minutes/cycle) were used in the 
simulations. Results are shown in Table 6, which confirm the claim. 
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Table 6. Tracking Error Due to Collector Vibrations 


First simulation 

vibration freouency » 
vibration amplitude « 
maximum steady-state error » 


.001 hertz 16.6 min/cycle) 
1 degree 

.0055 rad. at motor 
,002 arc-nin at antenna 


Second simulation 

vibration frequency = 
vibration amplitude = 
maximum steady-state error = 


.01 hertz {= 1.66 min/cycle) 
1 degree 

.55 rad. at motor 


.19 arc-min at antenna 


4. APPROACHES FOR MODELING THE FLEXIBLE COLLECTOR BODY 


It is recognized that the structural dynamics of the flexible collector 
body may be modeled by transforming it to an equivalent orthotropic plate. 
Then the wealth of well known theory of plates can be applied readily. 
However, a thorough study of the structural responses of the collector 
plate requires a major effort. The discussion here will focus on the 
two approaches v/hich may be used to obtain the desired model. They are the 
analytical approach and the simulation approach. 

The Analytical Approach In this approach the structure of the collector 
truss is first transformed to an equivalent grid, which is then trans- 
formed to an equivalent orthotropic plate as shown in Figure 32. The 
differential equation for an orthotropic plate, including the inertial- 
force tern, has the form [5] 



where u » u(x,y,t) is the deflection in the z-direction as shown in 
Figure 32(c), q is the external loading per unit area, h is the thickness, 
and 0 is the density of the plate material. The problem of modeling, simply 
stated, is the determination of constants D , D , and H. 

X y 

Transformation from grid work to an orthotropic plate has been 
thoroughly studied by Timoshenko and his cohort [1]. Let the geometry 
of a grid be defimed in Figure 33, then the constants for an equivalent 
orthotropic plate are given by 


D 


X 





t 

a 


2 
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(S2) 


where E is the modulus of elasticity, and and C 2 are torsional 
rigidities of the beams parallel to the x and y axes, respectively. When 
the beams have rectangular cross-section as shown, c^ and Cg are given by [6] 


_ G ,,3 ufl ^ 1 T 1 / nff h \ 

- 7 - t7» 


192 ^1 - 


nit h 


C2 * y 1^2 


192 ^2 


(83) 


“T TT i T 

11 ^ " n»l,3, n= ^ h 


5.— 


For most practical cases the infinite series converges rapidly; usually 
two or three terms are sufficient to give accurate result. For example, if 

* 1 , then 

^1 

y — ^ tanh 5— 

n-1,3,5 — h^ 

■ tanh Y '*■ "y tanh tanh . 

3 5 

« 0.91715 + ^ (0.9998) + jj^(D + - - - 


* 0.91715(1 + .004486 + .0003489 + - - - ) 


- 0.9216 


Similarly, the results of the series becane 1.00073 and 1.00414 for 

equals 2 and 4 respectively. Using the result for » 2 in (83) o'ves 
t^ ^1 
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Cl ■ j h[l - ij) (1.00073)] 

TI 

. I (^1)3 (0.6861) 

- I (j) (0.6861) 

« 0.0286 G 


Clearly C 2 *c^ if t 2 *t.j. The bending moments and twisting manents are 
related to the deflection u> by the following expressions. 


M 

X 


Eh^ ^1 

IT bT 
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2 

U) 



Ch3 ^2 
Eh 2 3 oi 



c 

b 


i ^ 

1 3x3y 


(84) 



c 

a 


2 ^ 

1 3x3y 


Analytic methods for transforming a truss structure as shown in Figure 
32(a) to an equivalent grid, or directly to an equivalent orthotropic 
plate have not been developed. It is this part of the study which requires 
a major effort. However, some related cases have been studied, such as 
plates with many holes [7-10] and corrugated pipes. [7-11]. The 
equivalent solid plate or smooth pipe will produce the same strain, defor- 
mation, or strain energy as the actual perforated plate or corrugated 
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pipe. The results for the actual structure have been carried out by testing, 
numerical calculation through point-matching or collocation technique, 
and other appropriate analysis such as energy method. 

The Simulation Approach Constant D , D , and H of Eq. (30) for a 

X y 

given truss nay be obtained by a pa’'ameter matching technique with the 
help of computer simulation. Well known computer programs exist, such 
as NASTRAN and STRDYN, which use finite element approximations for 
flexible structure. 


5. COMCLUSION 


A study of the pointing control system for SPS was accomplished. The 
study Included two major areas, namely, the analysis of a pointing error 
sensing method for the SPS antenna, and the design of a mechanical pointing 
control system for the antenna. In addition, approaches for modeling the 
flexible body collector dynamics were given. 

In the area of point error sensing, an analysis was made on a pre- 
viously proposed method for developing pointing angle signals to be used 
by a servomechanism to point the antenna toward a rectenna on ground. From 
the analysis result, the method was found to be feasible. A normalized 
technique was added to the method to produce error signals that are In- 
dependent of the received signal magnitude. Effects of vibration and noise 
on pointing control accuracy were analyzed and presented In the form of 
charts . 

In the area of mechanical pointing control of the antenna, a baseline 
satellite configuration was selected as shown in Figure 16 (p. 30) together 
with a set of baseline parameters. Analysis were made on various phenomena 
which might influence the accuracy of the pointing control. They include 
the nonlinear slip-ring friction of rotary joints, motor characteristics, 
effect of centrifugal force induced torque, effect of variation in antenna's 
moment of inertia, effect of finite collector mass on pointing control, and 
the effect of vibration modes of the collector. The friction and backlash 
characteristics of the gear train was not available for analysis and were 
ignored. A mathematical model of the plant, whose output is the antenna 
pointing direction, was developed. The pointing control system was designed 
using an optimal control approach with the help of the root-square-locus method. 
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The mechanical pointing accuracy, with the various anomalies present, was found 
to be well within 1 arc-minute. 

From the results of this study It appears that the desired pointing 
accuracy of ^ 3 min Is possible. Suppose a reasonable error budget Is 
+ 1 min for each of the major error sources of (1) noise, (2) vibration, 
and (3) mechanical effects (gear backlash, stricktion, etc.). Suppose 
further that we require the ♦ 3 min to correspond to at most SOS of the 
beamwidth 0g so that operation Is always In the linear part of the error 
pattern. Then from Figure 13 and 6g » 12 min we require M<3 subarrays In 
the linear array. The allowed vibration error Is therefore 1 m1n/12 min ■ 

0.033 or 8.3% of 9g maximin. From Figure 12 this allowance Is not exceeded 
for M«3 even for mode 2 (k»2) vibration (ai * 1 m). The noise error of 
+ 1 min can be achieved by providing sufficient power In the ground station 
pilot transmitter. 
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APPENDIX A. CSMP FOR THE FIRST ROTARY JOINT CONTROL SYSTEM. 


••••CONTINUOUS SYSTEM MCOELING 
•••PROliLEM input ST ATEMEMT S**^ 

• STEP RESPONSE WITH DRAG OISTURBANCE 
R«175. 

TO-IOO.O 

E-R-C 

X1«U.75*E 
X2*INTGRL(0.0yXl } 

X3«LEDLAC(3.164,0.193,X2) 

XA«DcRIV(0.0,X3) 

X5«2.36*X4 

V«X94X3 

E3-V-C2 

C2»A3.4*CI 

C1-OFRIV{0.0,C) 

C2»A3.A*C1 

X5»96.05*E3 

X6«T0^X5 

X7-0.00010173*X6 

X8olNTCRU0.0,X7> 

C-INTCRL ( 0.0,X8) 

• 

TIMER FINTlM-20.f OUTOEL-.2, OELT-.OOl, PR0EL-.2 
pRTPlT Ef V, C 

LA3EL STEP RESPONSE WITH DRAG DISTUP9ANCE 

END 

STOP 
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APPENDIX B. CSMP FOR THE SECOND ROTARY JOINT CONTROL SYSTEM. 


t^^^CC.'JTI'JUuUS SYSTEM niOOELING PROGRAM**** 

***PRC3LEV INPUT STATEMENTS*** 

* POINTING CCrNTROL WITH SHIFTING HYSTERESIS TYPE 

^ CF NONLINEAR FRICTION TORQUE 

* NF=-.573, DELT=.001, D=KF*DELT 
INITIAL 

COriSTANT SGNP = C. rTH= I 00 . , TDP=0. , TD*0 . , D=- . C00573 
DYNAMIC 
SORT 

R=350.*STEP( O.C)-2GO.*STEP (10.0) 

E=R-C 

Xl = 25.1t!*E 
X2=INTGRL(0.0,X1) 

X3=LEULAC(2.20,0.270,X2) 

XA-=DERIV{0.0,X3) 

X'5=2.36fX4 
V=X9+X3 
E3=V-C2 
C2=43.4*C1 
C1=0ERIV(0,U,C) 

C2=43.4*C1 
X5=96. 05*E3 
X6=TD+X5 

X7=0. 000 10173’«=X6 
RATE=INTGRL(0.0,X7J 
C=INTGRL (0.0, RATE) 

NOSORT 

SGN=SIGN{l.O,RAT£) 

IF (SON .£Q. SGNP) GOTO 30 i | 

TDP=0. I 

302 TO=TDP+D*R ATE I 

IF (A5S(TD) .GT. TM) TD=S I GN ( TM, TD) j 

GO TO 304 1 

301 IF (ABS(TDP) .EQ. TM) GOTO 303 i 

GOTO 302 ' 

303 TD=TOP 

304 TOP=TD 
SGNP=SGN 

SORT 

terminal 

timer FINTI ‘><=20. , 0L'T0EL=.2, DELT=.001, PRDEL = .2 
METHOD RKSFX 
PRTPLT C, E, RATE, TD 
LmBEL GTE? RESPONSE 
£ ND 
STOP 


